Maximal Quantum Violation of the CGLMP Inequality on Its Both Sides 
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Abstract 

We investigate the maximal violations for both sides of the d-dimensional CGLMP inequality 
by using the Bell operator method. It turns out that the maximal violations have a decelerating 
increase as the dimension increases and tend to a finite value at infinity. The numerical values are 
given out up to d — 1 6 for positively maximal violations and d = 2 x 10 5 for negatively maximal 
violations. Counterintuitively, the negatively maximal violations tend to be a little stronger than the 
positively maximal violations. Further we show the states corresponding to these maximal violations 
and compare them with the maximally entangled states by utilizing entangled degree defined by 
von Neumann entropy. It shows that their entangled degree tends to some nonmaximal value as the 
dimension increases. 
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Bell inequality [1], arising from the Einstein-Podolsky-Rosen (EPR) paradox [2], has developed into an effective 
criterion between the classical and quantum theory. The famous Clauser-Hornc-Shimony-Holt (CHSH) inequality 
[3] for two entangled spin-1/2 particles has always provided an excellent test-bed for experimental verification of 
quantum mechanics against the predictions of local realism. In 2002, two research teams independently developed 
Bell inequalities for high-dimensional systems: the first one is a Clauser-Horne type (probability) inequality for two 
qutrits [4, 5]; the second one is a CHSH type inequality for two qudits (d is an arbitrarily high dimension), now known 
as the Collins-Gisin-Linden-Masser- Popoescu (CGLMP) inequality [6]. The tightness of the CGLMP inequality was 
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demonstrated in Ref. [7]. Such an inequality is later modified into more compact forms [8, 9] and its numerical 
violations for high-dimensional situations have also been investigated in [10, 11]. 

The CGLMP inequality for two <i-dimcnsional systems is constrained by two classical bounds at its two sides, one 
being positive and the other negative [9]. Both two classical bounds can be violated by quantum behaviors. Such a 
property provides us the chance to inspect the quantum nonlocality in the high-dimensional cases. In principle, the 
maximally violating values of the inequality can be achieved by searching all possible measurements on all kinds of 
entangled states including mixed states. However, the complexity of numerical simulation will increases in power as the 
dimension increases, mainly because the number of variable parameters during the searching process is proportional 
to d 2 + d (for pure states). In Ref. [10], Chen et al. have investigated the maximal quantum violation of the CGLMP 
inequality at its positive side by using the Bell operator method. Such a method has a merit that it transforms the 
maximal violation problem into a pure mathematical calculation of solving the maximal and minimal eigenvalues 
of a matrix. By doing this, it needs a prerequisite that a set of phase rules for fixing the measurement projectors 
must be known so as to explicitly give out the matrix elements. In general, the determination of these rules needs 
a large quantity of numerical simulations. As it is, the rule for obtaining the positively maximal violation of the 
CGLMP inequality has been known before [6, 10], but the one for the negatively maximal violation is still unknown. 
Physically, it is necessary for us to examine whether the negative violation is larger than the positive side, since a 
stronger quantum nonlocality may be indicated while the positive violation would be enough if that is not the case. 

In this paper, we first simplify the Bell operator method appeared in [10] and then give out, respectively, the phase 
rules for obtaining the positively and negatively maximal violations of the CGLMP inequality. After that, the matrix 
elements of Bell operator are derived. Numerical results about maximal quantum violations are achieved by solving 
the maximal and minimal eigenvalues of the Bell operator. At last, the strength of quantum nonlocality is compared 
between the positive and negative violations, and their corresponding states are also discussed. 

The CGLMP inequality as an extended version of the Bell-CHSH inequality for c?-dimensional Hilbert space has 
the form [8] 



Id — Qn + Q12 — Q21 + Q22, 
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FIG. 1: The diagram of unbiased multiport beam splitter. In front of the ith input port of the device a phase shifter is placed 
to change the phase of the incoming photon by <p(i). The whole action can be described by unitary operator U(<j) a ) and U(f>b)- 



where Qij are the Bell-type correlation functions defined by probabilities in the way 
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with the spin S = (d—l)/2 and eigenvalues of the correlation function, / lJ '(m, n) = S— M(e(i— j)(m + n), d) [s(x) = 1 
for x > and — 1 for x < 0; M(x,d) = x mod d.]. In classical theory, the Eq. (1) has two bounds [9] as 



which can be violated in quantum theory. Following [4, 6, 10, 13, 14, 16], we use the unitary transformation setting 
of unbiased multiport beam splitters (UMBS) [12, 13] instead of a complete d-dimensional unitary transformation to 
define projectors [see Eq. (5)]. Unbiased d-port beam splitter (see Fig. 1) is a device with the following property: 
if a photon enters any of the d single input ports, its chances of exit are equally split among the d output ports. In 
fact one can always build the device with the distinguishing trait that the elements of its unitary transition matrix 
T are solely powers of the root of unity 7 = exp(i27r/d), namely, Tu — ^7 fc '- I n front of ith input port of the 
device a phase shifter is placed to change the phase of the incoming photon by <j)(i). These d phase shifts, denoted 
for convenience as a "vector" of phase shifts (f> — (0(0), 0(1), • • ■ , 4>{d— 1)), are macroscopic local parameters that can 
be changed by the observer. Therefore, unbiased <i-port beam splitter together with the d phase shifters perform the 
unitary transformation U((f>) with the entries Uki = Tki exp[i(f>(l)]. The diagram of such unitary operation has been 
displayed in the Fig. 1. 

An arbitrary entangled state of two qudits after Schmidt decomposition reads 
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with Schmidt coefficients {a*} being real and satisfying normalization condition \cti\ 2 = 1 [17]. One can see that if 
we take on = l/Vd for all i it is just the maximally entangled state. The quantum prediction of the joint probability 
P(A a = k,Bb = I) when A a and Bt are measured in the initial state |$) is given by 



P(A a = k,B b = l) 
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where life = |fc)(fc|, ft; = |Z)(i| are the projectors for observers A and B not including UMBS transformation, respec- 
tively. Substituting Eq. (5) into the high-dimensional Bell inequality (1), one gets the Bell expression for the state 
|*): 
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By using the summation formula 
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we can reduce Eq. (6) to a real- valued expression, i.e., 
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From Eq. (7), we can see that I4 depends on the state coefhcients {an} and phase shifts {0i,2(j), Vi,2(to)}. For 
the positively maximal values of Id, there have been several investigations: when |<j>) is restricted in maximally 
entangled states, the quantum violation was shown numerically in three dimension in [4] and in an arbitrary dimension 
in [6]; when |<j>) is an arbitrarily entangled state, the maximal quantum violations were displayed in [10, 11, 14]. 
For both cases, the latter one has a larger value of Id than the former and hence indicates a stronger quantum 
nonlocality. Particularly, the approach developed in Ref. [10, 14] is related to the Bell operator, which can be defined 
by reparametrizing 7^(1$)) according to 

7 d (|$))=Tr(B|$)($|) = ($|B|$), 

where B is the so-called Bell operator [15]. Under the bases {|jj)}, the non- vanishing elements of the Bell operator, 
Bj m = (jj\B\mm), from Eq. (7) are those off-diagonal elements with 
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If we know phases {0i,2(i), <£i,2(to)}, the whole matrix (Bj m ) then can be determined. Thus, the positively maximal 
or negatively minimal quantum violations of Eq. (3) would equal to the maximal or minimal eigenvalues of the matrix 
(Bj m ) and their corresponding eigenvectors are the states |$) that are capable to make the inequality maximal or 
minimal. As a result, a physical problem of searching the positively maximal or negatively minimal quantum violations 



of Eq. (3) becomes a mathematical problem of solving the maximal or minimal eigenvalues of the matrix (Bj m ) for 
all possible choices of phase shifts {<p a , <fib}- 

In order to determine the maximal violation at both hand sides of the Bell inequalities (1), we must first choose 
such an appropriate form of {4> ai ipb} as to Id can have maximal quantum violations. Just as pointed out in the Ref. 
[10, 14, 16], the positively maximal quantum values of Id can be achieved when one takes phase shifts by 



<t>i{j) = 0, = <pi(j) = ^, tp 2 (j) = -7^. (9) 



As to the negatively minimal values of Id, we perform the numerical search for the inequality and find that it can be 
achieved by piecewise taking 
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for < j < [^p] (square brackets here denote taking integer part value), and taking 

for [4=2] + l<j <d-2- [4=2] , and taking 

, / \ „ i / '\ (2d + j)n , .. (2d + j)ir . . (2d+7')7r 

0l(?)=O, fe(j) = / . <PlU) = 2 / . ^(j)=- V , (IOC) 

for d — 1 — [pp] < i < — 1- From Eqn. (9), (10a), (10b) and (10c), we can write the phase shifts into an uniform 
format 
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where rij take different values described above. The non-vanishing matrix elements of (Bj m ) under condition (11) 
can be expressed as 



7 

When rij = j, the Eq. (12) returns to Bj iTn — 2(1 — 5j m )/(d — 1) cos[(j — m)Tr/2d] appeared in [10]. Now that matrix 
elements of (Bj m ) are all known, the remainder is to get the extreme eigenvalues. As follows, we apply the so-called 
power method in the numerical computing field to achieve such task. At the same time the normalized eigenvectors 
corresponding to these extremum are just Schmidt coefficients {cti} of Eq. (4). 

It is natural to ask about the optimality of the chosen set of measurements. Certainly, the maximal and minimal 
values obtained by using these phase settings are only when the unitary transformation is restricted in the experimental 
settings of UMBS. For a general case, we have also performed a numerical search for the inequality, by using universal 
SU(d) transformation instead of the above UMBS settings up to d = 9. The maximal and minimal values of Id are in 
accordance with the results obtained by the above configuration and the numerical values have been shown in Fig. 2. 

In order to get a more absolute measure of nonlocality for high-dimensional Bell inequalities, let us consider the 
initial entangled state (4) mixed with some amount of noise [13] 



p(F) = Fp noise + (l-F)\$)($\, (13) 

where the positive parameter F < 1 determines the "noise fraction" within the full state and p n oiso = 1 /d 2 (the bold 
1 is a d x d unit matrix). The threshold minimal F m j n , for which the state p(F) still allows a local realistic model, will 
be our numerical value of the strength of violation of local realism by the quantum state | $) . It has an indication that 
the higher F m i n is, the more noise admixture will be required to hide the nonclassicality of the quantum prediction, 
and so the stronger the quantum nonlocality is. For both the right-hand side and the left-hand side of Eq. (3), we 
have the relations between F m j n and the maximal violations (/™ ax )$ for pure states 



^f n ht = i-2//r x a*», (14) 

//f n (l<f»- 



pleft _ i i 9 
min 1 



hd + (i - M 



d- 1 



As for high-dimensional Bell inequalities, the states |$) corresponding to maximal violations are generally not the 
maximal entangled states any more. In order to measure the deviation of |$) from maximal entangled states by their 
entangled degree, we use the von Neumann entropy of the reduced pure states (party A for example) pa — Ttb(|$)($|) 
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[17] which quantifies the bi-particle entanglement and is defined by 

d-l 

S(p A ) = —Tt(pa log pa) = -J2 a ^°& a l ( 15 ) 

i=0 

where the logarithms are taken to base two. For a pure state pab of two subsystems A and B, the von Neumann entropy 
of the reduced density operators Eq. (15) can effectively distinguish classical and quantum mechanical correlations, 
i.e., quantify entanglement [18]. For the maximal entangled pure state (or GHZ state) |V>g) = ^ Sj=o \33)i tne ^q- 
(15) for such state reads S ma x = S(pa) = — J2j=o(^/d) log(l/d) = log d that corresponds to the maximal entropy logd 
in a d-dimensional Hilbert space. So we can use the rate S(pA)/S max as a useful criterion to measure the deviation; 
the smaller the rate is, the bigger the deviation is. 

In Fig. 2 we show numerically the maximal violations up to d = 10 6 for positively maximal violations and d — 2 x 10 5 
for negatively maximal violations. It is seen from Fig. 2(a) (d) that the strength of quantum nonlocality denoted by 
F m i n have a decelerating increase with the growing system dimension d and it approaches a finite maximal value when 
d goes to infinity. The increasing trend also shows a stronger quantum nonlocality in a higher dimensional system. An 
interesting phenomenon is that although the quantum nonlocality becomes stronger in a higher dimensional system, 
the corresponding state, from Fig 2.(b)(e), has a decreasing entangled degree and hence deviates farther from the 
most entangled state. From Fig. 2(c) (f), these Schmidt coefficients of states corresponding those maximal quantum 
violations have a symmetric distribution and they obey the same pattern no matter how big the dimension is. By 
comparing Fig. 2(a-c) with Fig. 2(d-f), there are several differences between the positively and negatively maximal 
violations: (i) at the same dimension, the negatively maximal violation of Id has a little larger F m ; n than the positively 
maximal violation (e.g. when d = 2 x 10 5 , the negatively maximal violation has F m ; n = 0.487 while the positively 
maximal violation has F m ; n = 0.477); (ii) compared with the positively maximal side, there is a smaller deviation 
for the states corresponding to negatively maximal side (e.g. when d = 1 x 10 5 , the negatively maximal violation 
has S(pA)/S max = 0.848 while the positively maximal violation has S(pA)/S max = 0.826). Note that for d = 3, the 
negative side of Id has no violation with F m i n = and so at this point it does not satisfy the above general rules. 

In summary, continued with the previous paper [10], we have investigated the maximal violations of the CGLMP 
inequality for two d-dimensional systems further. The maximal violation occurs at the nonmaximally entangled state, 
which is the eigenvector of Bell operator with the maximal eigenvalue. Numerically, we display the maximal violating 
values up to d — 10 6 for the positively maximal violation and d = 2 x 10 5 for the negatively maximal violation. It 
gives us a straightforward picture on the trend of quantum nonlocality with the dimension of systems increasing. That 
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FIG. 2: (Color online)The violation illustration of the CGLMP inequality in the d-dimensional Hilbert space, (a), (b), (c) 
describe the negative maximal violations, while (d), (e), (f) describe the positive maximal violations. 
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is, the maximal violations increase as the dimension d grows and approaches a finite value when d becomes infinity. 
In addition, deviations of the entangled pure states that are in accordance with maximal violations from maximally 
entangled states have also been discussed. It is shown that as the dimension d increases, the deviation also approaches 
some constant limit. 

Particulary, by comparing the violations at both sides of Li, we show that the negatively maximal violation has a 
little stronger quantum nonlocality than the positively maximal violation has. This may be counterintuitive to the 
previous notion that the positively maximal violation of the CGLMP inequality is enough. Also, by giving out the 
phase rules and the elements of the Bell operator matrix, we have in fact given a convenient way to calculate the 
maximal violation of the CGLMP inequality on its both sides. It can serve as a useful reference for experimental setting 
when testing the CGLMP inequality. The experimental test has been performed to verify the CGLMP inequality 
for the first few high-dimensional systems [19], and it shows that indeed there exist nonmaximally entangled states 
violating the inequality more strongly than the maximal entangled ones . Thus, the nonmaximally entangled states 
in high-dimensional quantum systems as discussed in this paper may turn to be comparably useful when applied to 
quantum cryptography and quantum communication complexity [20, 21, 22, 23, 24], in which previous researches are 
mostly based on maximally entangled states. 
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Figure legends 

Figure 1: The diagram of unbiased multiport beam splitter. In front of the ith input port of the device a phase 
shifter is placed to change the phase of the incoming photon by 4>{i). The whole action can be described by unitary 
operator U(<p a ) and U(ipb)- 

Figure 2: (Color online)Thc violation illustration of the CGLMP inequality in the d-dimcnsional Hilbcrt space, 
(a), (b), (c) describe the negative maximal violations, while (d), (e), (f) describe the positive maximal violations. 
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